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Abstract 


Jet precession has previously been proposed to explain the apparently repeating features in the light curves of a few 
gamma-ray bursts (GRBs). In this paper, we further apply the precession model to a bright GRB 220408B by 
examining both its temporal and spectral consistency with the predictions of the model. As one of the recently 
confirmed GRBs observed by our GRID CubeSat mission, GRB 220408B is noteworthy as it exhibits three 
apparently similar emission episodes. Furthermore, the similarities are reinforced by their strong temporal 
correlations and similar features in terms of spectral evolution and spectral lags. Our analysis demonstrates that 
these features can be well explained by the modulated emission of a Fast-Rise-Exponential-Decay (FRED) shape 


light curve intrinsically produced by a precessing jet with a precession period of 1824507 


s, a nutation period of 


11.1*02 s and viewed off-axis. This study provides a straightforward explanation for the complex yet similar multi- 


episode GRB light curves. 
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1. Introduction 


Regardless of its different types of origin, which can be 
either the collapse of a massive star (Paczynski 1986; Woosley 
1993; Woosley & Bloom 2006) or the merger of binary 
compact stars (Eichler et al. 1989), a gamma-ray burst (GRB) 
central engine is believed to resemble the same accretion 
system which consists of a central object, an accretion disk, and 
a relativistic jet. In particular, if the central object is a black 
hole (BH), the angular momentum direction of the BH and the 
accretion disk can differ due to the anisotropic explosions of its 
progenitor star. As the outer part of the disk has a sufficiently 
larger angular momentum, it will maintain its direction and 
drive the BH and the inner part of the accretion disk to precess 
due to the Lense-Thirring torque and the viscosity of the disk 
(Lense & Thirring 1918; Bardeen & Petterson 1975). The jet 
launched from the inner region of the disk will follow the 


rotating black hole to precess (e.g., Reynoso et al. 2006; Liu 
et al. 2010; Lei et al. 2012). Such precession can naturally 
cause the change of observer angle (calculated with respect to 
the moving direction of the ejected material, see Section 3.1). 
In some cases, when the precession period is shorter than the 
burst duration, and the jet's opening angle is small enough, 
precession can affect the observed light curve by introducing 
periodic-like or missing emissions. Several previous attempts 
(e.g., Lei et al. 2007; Liu et al. 2010) have been made to 
correlate those features with observations. 

In light of previous studies, we sought to find additional 
GRBs that display those characteristics that may be attributed 
to precession. GRB 220408B, a recent burst co-detected by 
Fermi (Bissaldi et al. 2022), Konus-Wind (Lysenko et al. 
2022), Astro-Sat (Gopalakrishnan et al. 2022), and GRID (this 
work), quickly caught our attention due to its multiple similar 
temporal episodes. In this paper, we first performed a detailed 
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analysis of GRB 220408B from the perspective of its light 
curve properties and spectral evolution (Section 2). Motivated 
by the similarities of the three episodes in light curve profile, 
spectral evolution, and spectral lags, we proposed to use a 
precession-nutation model to explain the observed properties of 
GRB 220408B (Section 3). The summary and discussion are 
presented in Section 4. 


2. Observation and Data Analysis 
2.1. The Data 


GRB 220408B triggered the Gamma-ray Burst Monitor 
(GBM) aboard the NASA Fermi Gamma-ray Space Telescope 
(Meegan et al. 2009) at 07:28:04.65 universal Time on 8 Apr 
2022 (hereafter Ty) Gpm). It is also the third confirmed GRB 
observed by GRID (short for Gamma-Ray Integrated 
Detectors), a low-cost project led by students aiming to build 
an all-sky and full-time CubeSat network to monitor high- 
energy transient sources, including GRBs, in low Earth orbits 
(Wen et al. 2019, 2021). To date, GRID has collected several 
confirmed GRBs as well as dozens of GRB candidates, of 
which the first is GRB 210 121A (Wang et al. 2021). In this 
work, we mostly utilize the Fermi/GBM data in consideration 
of its wide spectral coverage and high temporal and spectral 
resolution. As a result of a large separation angle (7260?) 
between the pointing direction of the detector and the GRB 
location, the GRID data of GRB 220408B suffer from low 
signal-to-noise ratio! (S / N) and, therefore, are only displayed 
in the top panel of Figure 1 as an illustration. 

We retrieved the time-tagged event (TTE) data set of GRB 
220408B from the Fermi/GBM public data archive." Two 
sodium iodide (Nal) detectors, namely n6 and n7, with the 
smallest viewing angles with respect to the GRB source 
direction, were selected for our analysis. Additionally, the 
brightest bismuth germanium oxide (BGO) detector, bl, was 
also selected as it extends to a higher energy range. These data 
were then processed according to the standard procedures 
described in Zhang et al. (2011) and Yang et al. (2022) to 
investigate the bursts temporal and spectral properties, as 
detailed below. 


2.2. The Three-episode Light Curve 


We plot the GBM and GRID light curves together in Figure 
1, using the same bin size of 0.325 s and the same alignment 
time, To, at To.Ggw. While both light curves are extracted from 
the same 15-350 keV energy range, the S/N of the brightest 
peak of the GRID light curve is about one-fifth of that of the 
GBM light curve due to the former's large off-axis angle. 
Nevertheless, the majority of the significant peaks on both light 

S—B 


IS The signal-to-noise ratio is defined as S/N = PE where S is the count 


rate in signal region, B is the background rate, o(B) is the variance of B. 
17 https: / /heasarc.gsfc.nasa.gov /FTP/fermi /data/ gbm /daily / 
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Figure 1. Light curves (top three panels) and accumulated light curve (bottom 
panel) of GRB 220408B. The bin size is set to 0.325 s for all light curves. The 
top panel shows the light curve within the energy range from 15 to 350 keV by 
combining the data from all four GRID detector units. The second panel shows 
the light curve of Nal detectors n6, n7 and BGO detector b1 of Fermi/GBM in 
the energy range of 15—350 keV. The third panel shows the light curve of Nal 
detectors n6, n7 and BGO detector b1 of Fermi/GBM in the energy range of 
10-25 keV. The bottom panel shows the accumulated counts of the light curve 
in the third panel. The three episodes are marked in orange, green, and purple, 
respectively. The gray solid (dashed) lines are drawn at 0% (5%) and 100% 
(9596) of the accumulated counts. The vertical red dashed lines in the bottom 
two panels represent the To interval. 


curves coincide, which reinforce the usefulness of CubeSat 
detectors for GRB research, even in non-ideal observational 
conditions. 

GRB 220408B exhibits an overall Fast-Rise-Exponential- 
Decay (FRED; Kocevski et al. 2003) profile while retaining a 
complex substructure characterized by three apparently 
separated emission episodes. Following the method in Yang 
et al. (2020a) and Yang et al. (2020b), the burst duration in the 
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Figure 2. (a) Light curve profiles of the three episodes are aligned according to the middle peak positions within red vertical dashed lines. Matched peaks are marked 
with green vertical dashed lines. (b) The correlation between Episode II and Episode I is delineated in orange, and the correlation between Episode II and Episode III is 
delineated in purple. The shaded areas represent the corresponding 3c uncertainties. 


standard energy range of 15-350keV is calculated as 
T90,15-350 kev ~ 30 s (see also Bissaldi et al. 2022), counted 
from Ty + 1.5 s to To + 31.5 s. Such a Ton range, however, does 
not cover the third episode, which starts at around Tọ + 55 s. 

We noticed that the third emission episode becomes 
particularly significant in low energies, as shown in the third 
panel of Figures 1 and 3, indicating a strong spectral evolution 
across the three episodes. We were thus motivated to recalculate 
the burst's Too in a lower energy range between 10 and 25 keV 
to be 62.5s (see the bottom panel of Figure 1), which more 
accurately conveys the burst timescale and the central engine 
activities (Zhang et al. 2014). Such an energy range is also 
utilized in dividing the burst into three episodes with a visual aid 
of the pulse structures, as colored in the third panel of Figure 1. 

Interestingly, the three emission episodes display striking 
similarities with each other in terms of duration, pulse structure, 
and spectral evolution. A more comprehensive analysis of that 
focus will be conducted in the rest of this section. 


2.3. Similarity in Overall Temporal Profile 


The similarity of the light curve profiles is illustrated in 
Figure 2(a), where the light curves'® of the three episodes are 


15 Those light curves are extracted in the energy range of 10-100 keV to 
improve S/N and binned to 0.1s to increase the visibility of the detailed 
Structures. 


first aligned according to the middle peak positions (red vertical 
dashed lines). Interestingly, such an alignment automatically 
results in several other peaks matching along (green vertical 
dashed lines). 

To further quantify the similarity, we calculated the 
correlation coefficients between any pair of the three light 
curves in Figure 2(b). The Pearson correlation coefficient is 
0.68 with a p-value of 7.60 x 107°! between Episode I and 
Episode II, and is 0.42 with a p-value of 1.15 x 10? between 
Episode II and Episode III. The strong correlations among the 
three episodes indicate that they may have the same physical 
origin, which could account for their similar shapes. 


2.4. Similarity in Multi-wavelength Behaviors 


We then divided the energy range between 10 and 800 keV 
into eight bands and extracted multi-wavelength light curves 
from Fermi/GBM Nal detectors n6, n7 and BGO detector b1 
using the method described in Liu et al. (2022). As shown in 
Figure 3, the profiles of the multi-wavelength light curves, 
including their characteristics, such as the peak time and width 
of the pulses, clearly evolve in accordance with increasing 
energy. Such an evolution is commonly observed in GRBs and 
often measured as a spectral lag, which refers to the delay of 
the arrival time of gamma-ray photons in different energy 
bands (e.g., Norris et al. 2000; Yi et al. 2006). Both positive 
(i.e., higher-energy photons arrive earlier) and negative lags, as 
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Figure 3. The multi-wavelength light curves of GRB 220408B. The orange, 
green, and purple blocks mark Episode I, Episode II, and Episode III, 
respectively. 


well as the positive-to-negative lag transitions (e.g., Wei et al. 
2017; Du et al. 2021; Liu et al. 2022), have been observed in 
some GRBs. 

Following the method described in Zhang et al. (2012) and 
Liu et al. (2022), we calculated the energy-dependent lags for 
all three episodes using the multi-wavelength light curve pairs 
in Figure 3. The results are shown in Figure 4. Interestingly, a 
positive-to-negative transition feature, with a roughly con- 
sistent trend, is observed in the lag-E relations in all three 
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episodes, which further confirms the similarity of the three 
emission episodes. 


2.5. Similarity in Spectral Evolution 


We performed both time-integrated and time-resolved 
spectral analyses over the periods of Episodes I, II, and III, 
respectively. Our time-resolved spectral analysis seeks to track 
the spectral evolution in as much detail as possible. To do so, 
we divided the burst duration into 37 time-resolved slices (see 
Figure 5 and Table 1), each containing sufficient photon counts 
(i.e., 20 counts per spectral bin; Zhang et al. 2018) to ensure 
statistical validity. Within each slice, we extracted the count 
spectra of GRB 220408B from Fermi/GBM Nal detectors n6, 
n7 and BGO detector b1 following the procedures described in 
Zhang et al. (2011); Yang et al. (2020a, 2020b, 2022) and Zou 
et al. (2021). Corresponding background spectra are acquired 
by applying the baseline method (Yang et al. 2020a; Zou et al. 
2021) to the time interval from To — 98 s to Ty + 134 s for each 
energy channel. The response matrices of the detectors are 
generated using the GBM Response Generator. ? 

For each slice, as well as each episode, we performed a 
spectral fit using the Monte Carlo fitting tool MySpecFit (Yang 
et al. 2022). AII the spectra are fitted by a cutoff power-law 
(CPL) model formulated as (Yu et al. 2016) 


—(a+2)E 
P 


N(E)-—AE?e E , (1) 
where o, A, and E, are the photon index, normalization 
coefficient, and peak energy, respectively. 

The results of our spectral fitting are presented in Table 1. 
Based on the ratio of Profile Gaussian likelihood to the degree 
of freedom (PGSTAT/dof; Arnaud 1996) statistics, our results 
indicate that the CPL model can adequately fit the spectra of 
the three episodes and time-resolved slices. Using the best-fit 
parameters of the CPL model, we plot the spectral evolution of 
GRB 220408B in Figure 5, along with the total light curve 
summing up the GBM detectors n6, n7, and b1 between 10 and 
1000 keV. Both o and E, exhibit strong spectral evolution, 
roughly consistent with the tracking behaviors as observed in 
other GRBs (e.g., Lu et al. 2012). 

Further analysis of the spectral evolution is carried out on an 
episode-by-episode basis. In Figure 6, we plot a and E, as a 
function of the flux, F, in 10—10,000 keV of each time slice for 
the three episodes. One can see those observable pairs are 
strongly correlated and consistently follow the same tracks of 
logE = 0.46'003log F + 4.957912 and a = 0.5600: log F+ 
2040. respectively (see also Table 2 for the fitting 
parameters). Such consistency suggests that the spectral 
evolution patterns of the three episodes are similar, despite the 
fact that the global spectra of the burst undergo a hard-to-soft 


I? https: / /fermi.gsfc.nasa.gov /ssc/data analysis /rmfit/gbmrsp-2.0. 10.tar.bz2 
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Figure 4. Spectral lags between the lowest energy band (10-25 keV) and any higher energy bands of the three episodes (Episodes I to III from top to bottom). The 
zero lag is shown with a dashed gray line. The horizontal error bars represent the range of energy bands, and the vertical error bars indicate the 1c uncertainties. 


transition, which is typically attributed to GRB central engine 
characteristics (e.g., Zhang et al. 2018). 


3. Model and Fit 


The similarities of the three emission episodes, as well as the 
overall FRED shape profile, appear to point toward a uniform 
origin that produces the observed gamma-ray emissions in a 
repeatable manner. A natural explanation for such features is 
that the GRB jet may precess while propagating outward from 
the central engine (Portegies Zwart et al. 1999; Lei et al. 2007; 
Liu et al. 2010). In this section, we further test this hypothesis 
by quantitatively fitting the observed data with a precession 
model. In addition to the precession itself, our toy model also 
considers that the nutation (Portegies Zwart et al. 1999) of the 
jet can contribute to the substructure of the light curves. 


3.1. The Precession—Nutation Model 


Considering the precession and nutation of a GRB jet, the 
observer angle, 0, defined as the angle between the jet 
propagating direction and line of sight (LOS), varies as a 
function of time. A GRB can be significantly observed only 
when 0 is less than the jet's half-opening angle. The periodic 
change of 0 may cause the jet to sweep across the LOS 
intermittently, which, when taking into account the intrinsic 
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Figure 5. The spectral evolution of the CPL model. The horizontal error bars 
represent the time spans, and the vertical error bars indicate the 1c uncertainties 
of the best-fit parameters. 
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emission profiles together, can lead to complex shapes of GRB 
light curves, sometimes with repeating (Portegies Zwart et al. 
1999; Lei et al. 2007) and emission-missing (Wang et al. 2022) 
features. 

Our model is illustrated in Figure 7. A GRB jet with a half- 
opening angle, Oj, is propagating along its direction of £j. An 
observer resides within the 6,..-cone with an off-axis angle, 0, 
with respect to £;.,. The jet precesses with an angular velocity of 
Wpre along the z-axis while its rotating axis is nutating with an 
angular velocity of wyy. The x-y plane is set accordingly so the 
Cartesian coordinate is centered at the GRB central engine. We 
also assumed that the intrinsic emission from the jet is shaped 
as a FRED function (Kocevski et al. 2003), namely, 


rcd 

t r d r t rl p rr 
k(t = Fm t y 2 
(0 (=) d+r zu (2) 


where f is measured in the laboratory frame (the jet's local 
frame), tm is the time when the flux reaches the peak, Fm is the 
peak flux, r and d are the power-law exponents for the rise and 
decay, respectively. 

We then derived the observed flux of our model based on the 
above configuration. The direct effect brought by precession 
and nutation is the change of 0 as a function of time, which can 
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Figure 7. The schematic sketch of the precession-nutation jet model. 


be calculated as 
A(t) E arccos(Fjet (t) i fobs) 
= arccos[cos(A(1)) cos i + sin(A(t))sini cos(7(t))], 
(3) 


where A(f) is the angle between fj. (t) and z-axis, i is the angle 
between fp; and z-axis. 7/(f) is defined as 


n(t) = Pier (t) m Oops: (4) 


where Øjet and Pons represent the separation angles of directions 
of the jet (Fjet) and observer (fps) to the x-axis, respectively. 
According to the kinematical description of the angular 
evolution of the jet resulting from the precession and nutation 
(Portegies Zwart et al. 1999), A and Øje can be expressed as 


Wpre 


A(t) = Ao + 


tan Ao cos(wnuź), (5) 


nu 


w . 
Pier (t) = Piet,0 + Wpret + — SiN(Wput). (6) 
nu 
The increasing time intervals between the three episodes 
may be due to the slowing down of jet precession, which is a 
natural outcome of energy dissipation. We thus assumed a 
power-law decay of the precession angular velocity as 


t — to » 
Upre (t) = Wpre,0 t , (7) 


c 
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Table 2 
The Best-Fit Parameters of Linear Models for a—log F, log E,—log F, and log E,—a Correlations 
Line Model Episode I Episode II Whole Burst 
k b k b k b 

a = klogF + b 0.70220 3.48*021 0.677034 3.391133 0.56001 214508 
log E, = klogF + b 0.47002 4.9903) 0.55018 5.508 0.46 003 4.95*017 
log E, = ka + b 0.06022 2.33011 0.307033 2.307031 0.49* 009 2.56006 
where £ is the decay index, tọ is the offset time when the jet 3.2. The Fit 


begins to precess, fc is a characteristic timescale. 

We assumed a conical jet with a half-opening angle of 0j. 
and no moving material outside the cone. According to the 
derivation in Salafia et al. (2016), the observed GRB flux M(t, 0 
(A) at viewing angle 0(f) in the laboratory frame can be 
described by 


M (t, 0(t)) = Fo(t, 0 = 0) 


1, A(t) < 0%, 
TOO — 05) 
h-E, <0 < 4 
4— 4/2 01/3 
1 D jet 
—| > on > A(t) > Ger, 
2| 0 + 95r 
(8) 
where F,(t, 0 = 0) is the observed intrinsic flux when the LOS 
is centered on the jet axis, Ofa = jer — D T= (1 — 8257-2? ig 


the Lorentz factor, @ is the dimensionless radial velocity of 
the jet (i.e., 8 = Vie/ C, c is the speed of light, vi, is the speed 
of the jet), and D is the Doppler factor defined as 


1 . 
= Tasos 89r For GRBs, the Lorentz factor, I’, is 


typically a few hundred. In this work, we fix I' to be  — 300 in 
consideration that the bulk Lorentz factor does not significantly 
vary during the prompt emission phase. We also verified that 
our fitting result is not significantly affected by different values 
of I. 

Furthermore, one needs to convert the time in the laboratory 
frame to the observer frame by (e.g., Zhang 2018) 


1 — 8 cos0 (t) : 
Des 
where At is a parameter for adjusting the time offset of the 
model light curve. 
Finally, the observed flux can be calculated by substituting 
Equations (2)-(7) and Equation (9) to Equation (8), which can 
be written in form of 


M (tovs) = M (tovs, P), (10) 
where P represents the parameter set as P = {Xo, 
i, No» C"pre,0» Wnu» e to, fc, tm, T, d, Fn, At, Diets Dj in which 
No is defined as 1o = djet,.0 — ovs. 


fobs = At, (9) 


The next step is to fit our model (Equation (10)) to the 
observed light curve. With a fixed parameter of I' — 300, the 
free parameter set P of our model consists of the following 14 
items: 


1. The jet's half-opening angle Oje. According to Ryan et al. 
(2015), the maximum of jet is smaller than 0.5 radians. 
Thus the prior of Oje is set as a uniform distribution 
between 0 and 0.5 radians. 

2. The initial precession angle Ag. The prior of A is set as a 
uniform distribution between 0 and 7/2 radians. 

3. The observer's polar angle i. The prior of i is set as a 
uniform distribution between 0 and 7/2 radians. 

4. The initial phase rj. The prior of 7 is set as a uniform 
distribution between —7 and 7 radians. 

5. The initial precession angular velocity wpre,o. As there are 
at least two precession periods within the ~70 s duration 
of the burst, the prior of Wwpre,o is set as a uniform 
distribution between 0.18 and 1 rads~’. Such a range can 
account for the decrease of wpre- 

6. The nutation angular velocity w,,. The prior of wnu is set 
as a uniform distribution between 0.18 and 2.5 rad s !. 

7. The precession angular velocity decay index £. The prior 
of € is set as a uniform distribution between 0 and 0.5. 

8. The offset time fy of precession angular velocity decay. 
The prior of tọ is set as a uniform distribution between 
—50 and 10 s. 

9. The characteristic timescale fc of precession angular 
velocity decay. The prior of fc is set as a uniform 
distribution between 0 and 15 s. 

10. The peak time tm of the intrinsic FRED profile. The prior 
of t is set as a uniform distribution between 0 and 10 s. 

11. The rise power-law exponent r of the intrinsic FRED 
profile. The prior of r is set as a uniform distribution 
between 0.3 and 5. 

12. The decay power-law exponent d of the intrinsic FRED 
profile. The prior of d is set as a uniform distribution 
between 0.3 and 5. 

13. The peak flux Fm of the intrinsic FRED profile. The prior 
of Fm is set as a uniform distribution between 3000 and 
7000 cts s '. 
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14. The time offset At of the model light curve. The prior of 
At is set as a uniform distribution between —2 and 2 s. 


We then performed the fit using a self-developed Bayesian 
Monte-Carlo fitting package McEasyFit (Zhang et al. 2015), 
which is based on the widely used Multinest algorithm (Feroz 
& Hobson 2008; Feroz et al. 2009). This package can explore 
the complete parameter space efficiently to find the reliable 
best-fit parameters and determine their uncertainties realisti- 
cally by the converged Markov Chains. The log-likelihood 
function can be calculated as: 


n MN . 2 
th = In L(M|P) — D» S nn 2 , (11) 


i=1 Ci 


where fops,; and P(t55,;) represent the observed time and count 
rate of the i” data point of the light curve. F(t;5,;) is calculated 
by summing up the count rate of the GBM detectors n6, n7, and 
bl between 10 and 1000 keV at fobs, with a bin size of 0.25 s 
(blue curve in Figure 8(b)), M(fos,; P) is the model flux 
calculated at fobs, c; is the error of F(tops,;) estimated using the 
Poisson parameter confidence interval (Gehrels 1986). 
Furthermore, the model is required to reproduce two significant 
peaks in Episode III. Such a condition is guaranteed by forcing 
LL to be — oo whenever M (təs, P) < 3G during the time 
intervals of 61.35 < fobs S< 62.38 and 72.08 € toy, < 72.5 s, 
where G = 96.8 ctss ! is the variance of the observed 
background rate. 


3.3. The Result 


The best-fit model and the corner plot of the posterior 
probability distributions of the parameters are shown in Figure 
8. Our model successfully fits the observed light curve with 
PGSTAT/dof = 809.2/626. The best-fit parameters are listed 
in Table 3. Our results suggest that the precession-nutation 
model can well explain the main features of the observed light 
curve and point to an intrinsic FRED shape emission with 
r = 1S d = 209505 fn = 5.76191} s and Fm = 
4743.38 0514 ctss ! produced by a precessing jet with an 
initial precession period of 18.4 + 0.2 s and a nutation period of 
11.1 +0.2 s. Such an intrinsic shape was modulated to be a 
periodic-like and missing pattern as observed in GRB 
220408B. 

Based on the best-fit parameters in Table 3 and Equation (3), 
we can calculate that the observer angle 0 is always smaller 
than the jet's half-opening angle Oje, suggesting that the change 
of 0 does not dominate the change of the laboratory frame 
observed flux. The lower limit of the jet's half-opening angle 
can be derived as Gyct,lolim = Amax + imax = 1.05 x 107? rad. 
On the other hand, the precession-nutation effect modulates the 
shape of the observed light curve mainly through the 
conversion of photons' observed time between the laboratory 
frame and the observer frame (i.e., Equation (9)) rather than the 
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direct influence to the laboratory frame intrinsic light curve 
(i.e., Equation (8)). As a result of the conversion of arrival time 
between the two frames, the number of arriving photons is 
redistributed in the observer frame. At certain times, the arrival 
of photons is more concentrated, which results in the peak 
structures in the light curve. 


4. Summary and Discussion 


This paper proposes that the observed three-episode feature 
of GRB 220408B can be explained by a precessing jet. Based 
on the similarities between the three episodes in light curve 
profile, spectral evolution, and spectral lags, we concluded that 
they may have the same origin and may be the result of jet 
precession. A jet-precession model can be successfully used to 
fit the light curve of GRB 220408B, which assumes a FRED 
shape light curve that precesses and nutates with slowing 
precession angular velocity. Our fit suggests that the photon 
arrival time change in different frames resulting from the 
precession jet plays a prominent role in shaping the observed 
light curve when a GRB is observed off-axis. 

Our best-fit model is a precessing jet with an initial 
precession period of 18.4*02 s and a nutation period of 
11.1482 s, based on which we can further estimate the 
parameters of the central engine. The analytic expression of 
the precession period P can be expressed as 


7/13 . 
P = 2793 a}! aa sd 
M. Mo sS 


—30/13 
| q36/13 S, (12) 


(Q) 


where a, and M are respectively the dimensionless spin 
parameter (0 < ax < 1) and the mass of the BH, and M and a 
are respectively the accretion rate and the viscosity parameter 
of the disk. Assuming the gamma-ray radiation originates from 
the BH hyper-accretion processes, the isotropic gamma-ray 
luminosity and the mass accretion rate M can be related 
through 


M = dg e Las, (13) 


where ņ includes the beaming effect and efficiency of 
converting accretion material to gamma-ray radiation. The 
redshift of GRB 220408B remains uncertain, but we can make 
a rough estimation by assuming it adheres to the Amati relation 
for Type II GRBs. This relation is given as 
log £j; = a log Eiso — b, where the rest frame spectral energy 
peak is E,,= p.obs/ (1 +z), and the total isotropic energy 
Eio = 4nd?*fluence/(1 +z). The values for the best-fit 
parameters, a = 0.38 and b = 17.76, are obtained from a 
recent study (Sun et al. 2023). Using these parameters, we find 
that the best-fit redshift of GRB 220408B corresponds to a 
value of z = 2.0. Then, we can calculate the corresponding 
isotropic gamma-ray luminosity L,;,, to be approximately 
4.5 x lO ergs '. By substituting Equation (12) into 
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i = 3.085922! x 10-4 


a 

— 

HI 

10730517888] c 

=] 

o 

O 

funu = 0.56753 v 

i 

oO 

c 

[^ 
pre,o = 0.342109] 
° À 
gor 
ty 7027919355 


SL 


— observed light curve 
— model 


Time (s) 


MI 


d. 
e 
D 


eee m» | € 


‘= pde» e e» je | «ei a> 


SFOS Z 


Sm] ID [DS 1S | dd €] aR 6) oh 


du elas de 


ry ET 

@ db s. 

we e^. 

$ ap wee$e^ 

E: Ni E: BS ISP ISS N o 
LE Te E 331: E a 


Figure 8. (a) Corner plot of the posterior probability distributions of the parameters. The red error bars represent the lc uncertainties. The contours in the 2D 
histograms represent lo, 1.50, and 2c uncertainties, respectively. (b) Light curves of GRB 220408B and the best-fit precession-nutation model. 


Best-fit 
0.33033 
6.45103 x 1074 
3.09002 x 1074 
3.05 50; 
0.347091 
0.570005 
0.327001 
0.28788 
11.267031 
5:16:05 
1.597904 
222! 009 
4743.381874 


Table 3 
The Best-Fit Parameters of the Light Curve in 10-1000 keV Range 

Parameter Range 
Üje. (rad) (0, 0.5] 
Ac(rad) [0, 7/2] 
i(rad) [0, 7/2] 
no(tad) [-7. 7) 
pre, o (rad s7!) [0.18, 1.0] 
Way(tad s7!) [0.18, 2.5] 
€ [0, 0.5] 
to (s) [—50, 10] 
tc (s) (0, 15] 
tm (S) (0, 10] 
r [0.3, 5.0] 
d [0.3, 5.0] 
Fn (cts s7!) [3000, 7000] 
At (s) [—2.0, 2.0] 


— 1.267385 


10 


Equation (13), we obtain 


13/7 —17/7 —30/7 
see (atc) es)" Gs 
18.4s 0.9 0.05 


. 30/7 —36/7 
x Lyis = . — (14) 
4.5 x 10°'erg s^! 0.0125 


Here, we utilized a ~ tpe ^ = = 0.01 (e.g., Hou et al. 
2014b) for a typical accretion timescale of t,.. = 80s, and a 
preferred a, = 0.9 for the collapsar model (e.g., Popham et al. 
1999; Lei et al. 2007; Liu et al. 2010). Additionally, we 
adopted a value of 0.05 for the 7 parameter, although it is 
essential to acknowledge that this choice is somewhat arbitrary, 
given the inherent uncertainty associated with this parameter in 
theory. In the literature, a range of values between 0.01 and 
0.15 have been employed (e.g., Cannizzo & Gehrels 2009; Hou 
et al. 2014a, 2014b). With the parameters mentioned above, the 
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estimated black hole mass of 2.1 Mọ aligns with predictions 
from collapsar models (Popham et al. 1999). 

Additionally, it is worth noting that our model maintains a 
fixed jet bulk Lorentz factor at the typical value of I' = 300. 
Our analysis has verified that variations in the I' value have 
negligible impact on the profile of the model's light curve, 
thereby not substantially influencing our fitting outcomes. 
Furthermore, in view of computation costs, our model, which 
already has 14 free parameters, does not incorporate the 
reproduction of the observed spectral evolution. Although, in 
principle, such an evolution can be attributed to the Doppler 
factor change predicted by our model, a realistic model should 
also take into account the intrinsic central engine behaviors that 
result in the observed spectral evolution, which adds even more 
complexity to the model but will be left for future study. 
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